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Within the Skyrme-Hartree-Fock (SHF) approach, we show that for a fixed mass number A, both 
the symmetry energy coefficient &sym (A) in the semi-empirical mass formula and the nuclear matter 
symmetry energy E synl (pA) at a subsaturation reference density pA can be determined essentially 
by the symmetry energy E Byul (po) and its density slope L at saturation density po. Meanwhile, we 
find the dependence of a sym (^4) on E syul (po) or L is approximately linear and is very similar to the 
corresponding linear dependence displayed by E synl (pA), providing an explanation for the relation 
E synl (pA) ~ a sym (A). Our results indicate that a value of E syln (pA) leads to a linear correlation 
between E synl (po) and L and thus can put important constraints on -E S ym(po) and L. Particularly, 
the values of E synl (po) = 30.5 ± 3 MeV and L = 52.5 ± 20 MeV are simultaneously obtained by 
combining the constraints from recently extracted E Bynl (pA — 0.1 fm -3 ) with those from recent 
analyses of neutron skin thickness of Sn isotopes in the same SHF approach. 
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I. INTRODUCTION 



The study of the nuclear matter symmetry energy 
E sym (p), which essentially characterizes the isospin de- 
pendent part of the equation of state (EOS) of asym- 
metric nuclear matter, is currently an exciting topic of 
research in nuclear physics. Knowledge about the sym- 
metry energy is essential in understanding many aspects 
of nuclear physics and astrophysics [ij-ij as well as some 
interesting issues regar ding possible new physics beyond 
the standard model [7H10|. In recent years, significant 
progress has been made in determining the density de- 
pendence of E sym (p) 0,0], especially its value E sym (po) 
and its density slope L at saturation density pq. While 
constraints on -E S ym(po) and L from different experimen- 
tal data or methods become consistently convergent [TT1 -- 
[T3 |. they are still far from an accuracy required for un- 
derstanding enough precisely many important properties 
of neutron stars [3, LL5| . To narrow the uncertainty of the 
constrains on E synl (po) and L by using more accurate 
data or new methods is thus of crucial importance. 

Recently, based on the calculations within the droplet 
model and mean field models using a number of different 
parameter sets, Centelles et a! found that the symme- 
try energy coefficient a sym (A) of finite nuclei with mass 
number A in the semi-empirical mass formula can ap- 
proximately equal to nuclear matter symmetry energy 
E &ym {pA) at a reference density pa in the subnormal 



density region, i.e., E sym (p A ) 



,(A) [16(. This re- 



lation provides the possibility to directly determine the 
symmetry energy at subnormal densities from the semi- 
empirical mass formula and also has many important 
implications for extracting the symmetry e nerg y from 
isospin dependent observables of finite nuclei [16| ■ While 
this relation has been used to extract information on the 
symmetry energy around the normal density fl6l [l7j , its 
microscopic explanation is still missing. 
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Within the Skyrme-Hartree-Fock (SHF) energy density 
functional, we demonstrate in the present work that both 
(A) and E sym (p A ) can be determined essentially by 
(po) and L and meanwhile they display very similar 
linear dependence on E syia (po) or L, and thus provid- 
ing an explanation for the relation E syul (pA) ~ a syr[L (A). 
Furthermore, we show that a value of E sym (pA) can put 
important constraints on £ , S ym(po) and L. Combining the 
constraints from recently extracted E sym (p A = 0.1 fm -3 ) 
with those from recent analyses of existing data on neu- 
tron skin thickness of Sn isotopes [l8j within the same 
SHF approach leads to stringent constraints simultane- 
ously on E sym (po) and L. 



II. 



E sym (p) AND a sym {A) IN THE 



SKYRME-HARTREE-FOCK APPROACH 

The EOS of isospin asymmetric nuclear matter, given 
by its binding energy per nucleon, can be expanded to 
2nd-order in isospin asymmetry S as 



E(p,6) = E {p) + E sym (p)6 2 + 0(<5 4 ) 



(1) 



where p = p„ + p v is the baryon density with p n and 
p p denoting the neutron and proton densities, respec- 
tively; S = (p n —Pp)/p is the isospin asymmetry; Eq(p) = 
E(p, S = 0) is the binding energy per nucleon in symmet- 
ric nuclear matter, and the nuclear symmetry energy is 
expressed as 



E sym (p) 



1 d 2 E(p, 5) , 
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with x = E r £s -. 

3,90 

and K sym = 9pg 



The coefficients L = 3po 



dE syla (p) | 



dp 1 



dp \P=Pa 

\p=p are the slope and cur- 
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vature parameters of the symmetry energy, respectively. 
Within the standard SHF approach, the symmetry en- 
ergy can be written as (see, e.g., Ref. fl9l ]) 

h 2 /3tt 2 \ 2/3 * I , 
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with 9 sym = 3tiXi - i 2 (4 + 5x 2 ) and a, t - t 3 . x — x 3 
being the Skyrme interaction parameters. 

As shown in Refs. [H|,[2(|, the 9 Skyrme interaction pa- 
rameters, i.e., (j, t — i 3 , xq — x 3 can be expressed analyti- 
cally in terms of 9 macroscopic quantities po, Ea(po), the 
incompressibility Kq, the isoscalar effective mass to* , 
the isovector effective mass m* , E sym (po), L, gradient 
coefficient Gs, and symmetry-gradient coefficient Gy. In 
terms of these macroscopic quantities, the symmetry en- 
ergy can be rewritten as 



with 



E sym (p) = AxEsymipo) + B X L + d, 



A 1 = (tu - «r)/( 7 - 1) 
B 1 = (uT-u)/[3( 7 -l)] 
C X = E^Jp Q )u^ + Du^ 

(37-2)^ y n m (p ) + (3 7 -5)£> 
3(7-1) 



(4) 



(5) 
(G) 
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where u = p/po is the reduced density; £*y™(po) = 
2 / 2 \ 2/3 

( ^i - / 90 ) * s ^ ne kinetic symmetry energy at po; and 
the parameters D and 7 are defined as [2f|, [22| 
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(9) 
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with G = ^< in and £° n = (|- 2 Po) 2/3 . 

The symmetry energy coefficient a sym (A) of finite nu- 
clei in the semi-empirical mass formula can be expressed 
as O 



fA , E sym (p Q ) , 9E sym (p ) 1/3 

a sy m(A) = — i - ; 1 Wltn ^ = 77^ ^ ' , 

1 + XA 4Q 

(10) 



where the Q parameter is the so-called neutron-skin stiff- 
ness coefficient in the droplet model [H, [H| and it is 
related to the nuclear surface symmetry energy pBl [26| . 
Usually for a given nuclear interaction, the Q parameter 
can be obtained from asymmetric semi-infinite nuclear 
matter calculations [ID, [2J, [26l - [28j . As a goo d ap proxi- 
mation, the Q parameter can be expressed as [27] 
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(11) 



where r nm = (^irpo) i/u is the radius constant of nu- 
clear matter and a is the diffuseness parameter in the 
Fermi-like function from the parametrization of nuclear 
surface profile of symmetric semi-infinite nuclear matter. 
Many calculations [25l427l | have indicated a m 0.55 fm 
and then 2a/r nm w 1. Therefore, the xa parameter can 
be approximated by 



xa = (L - A sym /I2) 



E sym {po) 



(12) 



It should be noted that Eq. (|f 2|) is a good approxima- 
tion for evaluating the a sym (A), and within the standard 
SHF energy density functional the difference between the 
value of a sym (A = 208) from Eq. (fT2"]) and that of using 
the exact Q parameter obtained from asymmetric semi- 
infinite nuclear matter calculations is essentially less than 
I MeV [2^]. Furthermore, within the standard SHF ap- 
proach, K xvr n ca n be written in terms of the macroscopic 
quantites as [22| 



K sym = 3 7 L + ^(p )(37-2) 

+2,0(5 - 3 7 ) - 9jE aym (p ), 



(13) 



and thus we have 
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For < f , a sym (A) in Eq. ([TO]) can be expanded as 
a sym (A) = E sym (p )(l -x A + x\ ). (15) 



Neglectiing the x A and higher-order terms in Eq. ([f5 
leads to 



(A) = A 2 E syni (p Q ) + B 2 L + G 2 , (16) 



with 
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However, the convergence of the expansion in Eq. (jT5|) 
is usually very slow and thus Eq. (fTB|) is a very bad 
approximation to a sym (A) even for heavy nuclei pa. I30I] . 
A much better approximation could be obtained by the 
two- variable Taylor expansion with respect to E sym (po) 
and L at a point of E sym (po) = Sq and L = Lq as 

a S ym(^l) = a sym(A)\ Eaym (p )=So,L=L 



+ {E sym (p ) - So) 
+ (L-Lo 



dE sym (p 

da syT a(A) 



and keeping only the first-order terms leads to 
(A) = A 3 E sym ( Po ) + B 3 L + C 3 , 



with 



and 



A 3 = 
B 3 = 
C 3 = 
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III. 



NUMERICAL RESULTS AND 
DISCUSSIONS 



One can see from Eq. (j4j that E sym (p) is linear func- 
tions of E sym (po) and L with the coefficients A\, B\ 
and C\ determined by the density p and nuclear matter 
macroscopic quantities po, Eo(po), Kq, to* , and to* . 
Meanwhile, from Eqs. (ITU1) and (HU), one can see that 
asym(^4) is determined by the mass number A and also 
the nuclear matter macroscopic quantities. In particular, 
a sym (A) can also be linear functions of E sym (p ) and L 
if the approximation (|16p or (|2"Tj) is valid. As mentioned 
previously, the relation a sym (^4) s=s £^ sym (/9^i) has been ob- 
served within mean field models using a number of differ- 
ent parameter sets for the nuclear effective interactions. 
Particularly, one finds a. 
fm~ 3 ), a sym (A = 116) ' 



' sym 



(A = 40) « £ sym (p 



sym (A = 208) w S sym ( PA 
« S syra (p A - 0.093 fin" 3 ), 
= 0.08 fm" 3 ) [111. This fea 



0.1 

and 



ture implies that a sym (A) and E sym {p A ) would display 
similar correlation with each nuclear matter macroscopic 
quantity among L, G v , G s , E (p ), E sym (p Q ), K , 
to* , and po , which completely determine the 9 Skyrme 
interaction parameters a, — t 3 , xq — x 3 . In the follow- 
ing, we show that this is indeed the case by analyzing the 
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FIG. 1: (Color online) E aym (p A ) with p A = 0.09,0.10, and 
0.11 fm -3 as well as a syul {A = 208) and its approximation 
with Taylor expansion in Eq. \21\ (with So = 30 MeV and 
L = 60 MeV) from SHF with MSL0 by varying individually 
L (a), G v (b), G s (c), E (po) (d), E Bym (p ) (e), K (f), m* j0 
(g), m* fi (h), po (i), and W (j). 



correlations of E sym (p = 0.1 fm -3 ) and a sym (A = 208) 
with the nuclear matter macroscopic quantities. We have 
also checked the cases of A = 116 and 40, and obtained 
the similar conclusion as in the case of A = 208 and con- 
firmed the relations a syin (A = 116) «J E sym (p A = 0.093 
fm- 3 ) and a sym (A = 40) « E sym (p A = 0.08 fm" 3 ). 

As a reference for the correlation analyses based on 
the standard SHF energy density functional, we use in 
the present work the MSL0 parameter set [18[, which is 
obtained by using the following empirical values for the 
macroscopic quantities: po = 0.16 fm , Eq(pq) = —16 
MeV, K = 230 MeV, to* = 0.8to, to* n = 0.7to 



"v,q 

E 8ym (p ) = 30 MeV, and L = 60 MeV, G v = 5 MeV-fm 5 , 
and Gs = 132 MeV-fm 5 . And the spin-orbit coupling 
constant Wq = 133.3 MeV -fm 5 is used to fit the neutron 
P1/2 —P3/2 splitting in 16 0. It has been shown [l8[ that 
the MSL0 interaction can describe reasonably the binding 
energies and charge rms radii for a number of closed-shell 
or semi-closcd-shcll nuclei. It should be pointed out that 
the MSL0 is only used here as a reference for the correla- 
tion analyses. Using other Skyrme interactions obtained 
from fitting measured binding energies and charge rms 
radii of finite nuclei as in usual Skyrme parametrization 
will not change our conclusion. 

Shown in Fig. [1] are E syia (pA) with p A = 0.09,0.10, 
and 0.11 fm™ 3 as well as a sym (^4 = 208) (from Eq. (pTJ|) 
with approximation in Eq. (|12p ) and its approximation 
with Taylor expansion in Eq. ([51]) (with So = 30 MeV 
and L Q = 60 MeV) obtained from SHF with MSL0 by 
varying individually L, G v , G s , E (p ), E sym (p ), K , 
to* , to* , poj and Wo, namely, varying one quantity at 
a time while keeping all others at their default values in 
MSL0. We note here that using the exact Q parame- 
ter obtained from asymmetric semi-infinite nuclear mat- 
ter calculations to evaluate the a sym (^4) does not change 



4 



our conclusions [29|. It is interesting to sec that, within 
the uncertain ranges considered here for the macroscopic 
quantities, a sym (A — 208) displays strong correlations 
with both E sym (po) and L while it is almost no de- 
pendence on other macroscopic quantities (From Eqs. 
(fTOll and Ifil). 



^(A) is independent of Gy, Gs, and 
Wo)- This is understandable since a sym (^4) is determined 
uniquely by the three lowest-order characteristic parame- 
ters of the symmetry energy, i.e., E sym (po), L, and K syB1 
as seen in Eqs. (fTUj) and (|12p while K sym has been found 
to strongly correlate with -E S ym(Po) and L but exhibit 
very weak dependence on other macroscopic quantities 
within the standard SHF energy density functional as 
shown in Ref. [2(|. Furthermore, a sym (A = 208) displays 
approximately linear correlations with both E sym (po) 
and L which is demonstrated by the good approxima- 
tion of Eq. (J5TJ to a sym (A = 208) observed in Fig. Q] 

Similarly, one can see from Fig. [Tjthat E syra (pA) with 
Pa = 0.09, 0.10, and 0.11 fin -3 display strong linear cor- 
relations with both E sym (po) and L while they are almost 
independent of other macroscopic quantities except with 
a small dependence on po (Note E sym (pA) is indepen- 
dent of Gy, Gs, and Wo). These results indicate that 
dsym(^4 = 208) and E sym (pA) at subsaturation reference 
densities pa — 0.09,0.10, and 0.11 fur 3 are essentially 
determined by E sym (po) and L. Furthermore, Fig. Q] 
shows that both a sym (A = 208) and E sym (pA) display 
very similar linear dependence on -E S ym(po) or L. Espe- 
cially, it is seen from Fig. [T]that E sym (pA) with pa = 0.10 
fm -3 gives the best fit to a sym (A = 208). In other words, 
for any Skyrme force determined by the 9 parameters a, 
to — t^, xo — £3 or equivalently the 9 macroscopic quan- 
tities L, Gy, G s , E (p ), E sym (p ), K , m* , m* j0 , po, 
the value of E sym (pA — 0.10 fur 3 ) will approximately 
equal to that of a sym (A = 208), and this explains the 
relation E sym (pA = 0.10 fur 3 ) w a sym (A = 208) ob- 
served within mean field models using a number of dif- 
ferent parameter sets. Furthermore, one can see that the 
possible small deviations between E sym (pA = 0.10 fm -3 ) 
and a svm (A = 208) observed for some different parame- 
ter sets [16j may be mainly due to the different values of 
L, E sym (p ), and/or p . 

Since E sym (pA = 0.10 fm -3 ) is essentially determined 
by E syirL {po) and L and displays linear correlations with 
the latters, a determination of E sym (pA = 0.10 fm -3 ) 
will then put important constraints on E syul {po) and L. 
As a matter of fact, the value of the symmetry energy 
around 0.1 fm -3 has been heavily under investigation in 
recent years in the literature |3ll446| . For example, an 
analysis of the giant dipole resonance (GDR) of 208 Pb 
with Skyrme forces suggests a constraint E sym (pA = 
0.10 fin- 3 ) = 23.3 - 24.9 MeV 01 while a relativistic 
mean- filed model analysis of the GDR of 132 Sn leads to 
E Bym (p A = 0.10 far 3 ) = 21.2 - 22.5 MeV Q. In a 
recent work [13, Liu et al. extracted the symmetry en- 
ergy coefficients a sym (A) for nuclei with mass number 
A = 20 — 250 from more than 2000 measured nuclear 
masses and they obtained a value of 20.22 — 24.74 MeV 
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FIG. 2: (Color online) Contour curves in the E synl (po)-L 
plane for E syln (p = 0.1 far 3 ). The region between the two 
thick solid lines represents the constraint obtained in the 
present work with 20.22 MeV < E sym (p A = 0.10 fm -3 ) < 
24.74 MeV while the region between the two thick dashed 
lines is the constraint from the SHF analysis of neutron skin 
data of Sn isotopes within a 2a uncertainty [l|| . The shaded 
region represents the overlap of the two constraints. 



for a syul (A = 208) within a 2a uncertainty and thus we 
have E sym (p A = 0.10 fm~ 3 ) = 20.22 - 24.74 MeV ac- 
cording to E syia (pA = 0.10 fm- 3 ) « a sym (A = 208). 
In the following, as a conservative estimate, we will use 
E sym (p A = 0.10 fm" 3 ) = 20.22 - 24.74 MeV since it is 
essentially consistent with the other two constraints from 
the GDR of 208 Pb and 132 Sn. 

Shown in Fig. [2]arc the contour curves in the E sym (po)- 
L plane for E sym (p = 0.1 fm~ 3 ) obtained from SHF with 
MSL0 by varying individually E sym (po) and L. The re- 
gion between the two thick solid lines in Fig. [^represents 
the constraint of 20.22 MeV < E sym (p = 0.10 fur 3 ) 
< 24.74 MeV from the nuclear mass restrictions on 
<^sym(A = 208). Also included in Fig. [5] is a recent con- 
straint (the region between the two thick dashed lines) 
from SHF analysis of neutron skin data of Sn isotopes 
within a 2a uncertainty (l8l |. It is interesting to note 
that the constraint from the neutron skin data of Sn 
isotopes suggests L decreases with increasing E sym (p ) 
while the constraint from E sym (p = 0.10 fm -3 ) displays 
opposite behaviors. This feature allows us to extract si- 
multaneously the values of both E sym (pQ) and L with 
higher accuracy, namely, E sym (po) = 30.5 ± 3 MeV and 
L = 52.5 ± 20 MeV by combining the two constraints as 
illustrated by the overlap of the two constraints repre- 
sented by the shaded region in Fig. [5J Furthermore, it 
is seen from Fig. [5] that a value of E synl (p = 0.10 fnr 3 ) 
puts a strong linear correlation between E sym (po) and L, 
i.e., L increases linearly with E syul (po), which has been 
extensively observed in the parameterization for nuclear 
effective interactions in mean-field models [HUMS- 
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The simultaneously extracted values of E sym (p Q ) = 
30.5 ± 3 MeV and L = 52.5 ± 20 MeV from the same 
SHF approach within a 2a uncertainty are essentially 
overlapped with other constraints extracted from dif- 
ferent experimental data or methods in the literature 
[II [ll Hill, [II 111 111 EE&Bi (seeRef. for a re- 
cent summary) but with higher precision. In particular, 
these extracted values are in remarkably good agreement 
with the E sym (p ) = 31.3 ±4.5 MeV and L = 52.7 ±22.5 
MeV extracted most recently from global nucleon optical 
potentials constrained by world data on nucleon-nucleus 
and (p,n) charge- exchange reactions [H[ ■ The extracted 
value of L = 52.5 ± 20 MeV also agrees well with the 
value of L = 58 ± 18 MeV obtained from combining the 
constraint from the neutron skin data of Sn isotopes fl~8T ] 
with those from recent analyses of isospin diffusion and 
the double neutron/proton ratio in heavy- ion collisions 
at intermediate energies fill ] . Furthermore, the extracted 
value of L = 52. 5 ±20 MeV is consistent with the value of 
L = 66.5 MeV obtained from a recent systematic analy- 
sis of the density dependence of nuclear symmetry energy 
within the microscopic Brucckncr-Hartrcc-Fock approach 
using the realistic Argonnc V18 nucleon-nucleon poten- 
tial plus a phenomenological three-body force of Urbana 
type [5(|. It should be stressed that the simultaneously 
extracted values of E sym (po) and L in the present work 
are obtained from the same Skyrme-Hartree-Fock energy 
density functional. 

IV. SUMMARY 

We have analyzed the correlations of the nuclear mat- 
ter symmetry energy E sym (p) at a subsaturation refer- 
ence density pa and the symmetry energy coefficient 
flsym(^l) of finite nuclei in the semi-empirical mass for- 
mula with nuclear matter macroscopic quantities within 
the Skyrme-Hartree-Fock energy density functional. We 
have shown that E sym (pA) displays explicitly linear cor- 
relations with E sym (po) and L and it is essentially de- 
termined by the latters but almost no dependence on 
other macroscopic quantities except a small dependence 
on the saturation density po- These features imply that 
a fixed value of E sym (pA) will lead to strong linear cor- 



relation between E syr[L (po) an d L. Furthermore, we have 
found that the two macroscopic quantities E sym (po) and 
L essentially determine the value of a syin (A) and the lat- 
ter displays approximately linear correlations with both 
E sym (po) and L. In particular, the correlation between 
flsym(^l) and E sym (po) (L) is found to be very similar to 
that between E syirL {pA) and E sym (po) (X), thus providing 
an explanation for the relation E sym (pA) ~ a sym (A) ob- 
served in mean field models using a number of different 
parameter sets. 

Using the relation E sym (pA) ~ a sym (A), we have 
demonstrated that within the Skyrme-Hartree-Fock 
energy density functional, the value of E sym (p = 
0.10 fm" 3 ) = 20.22 - 24.74 MeV extracted recently from 
nuclear masses within a 2a uncertainty [17| can put im- 
portant constraints on E sym (po) and L. Combining these 
constraints with those from recent analyses of existing 
data on neutron skin thickness of Sn isotopes based on 
the same Skyrme-Hartree-Fock approach within a 2a un- 
certainty [18[ allows us to extract simultaneously the val- 
ues of both E sym (p ) and L, i.e., E sym (p ) = 30.5±3 MeV 
and L = 52.5 ± 20 MeV. These extracted values are es- 
sentially consistent with other constraints extracted from 
different experimental data in the literature but with 
higher precision. 

In the present work, all analyses are based on the stan- 
dard Skyrmc-Hartrcc-Fock energy density functional. It 
will be interesting to see how our results change if dif- 
ferent energy-density functionals are used. On the other 
hand, it will be also interesting to sec how our results, 
especially the new constraints on E sym (po) and L ob- 
tained in the present work, can give implications for the 
neutron-skin thickness of heavy nuclei, the isovector gi- 
ant dipole resonance of finite nuclei, and properties of 
neutron stars. These studies are in progress. 
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